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I. INTRODUCTION 



Bound states of heavy quarks provide a powerful test- 
ing ground for quantum chromodynamics (QCD). Exper- 
imental studies of charmonium and bottomonium spec- 
troscopy have helped uncover some of the key aspects of 
the quarkonium potential |l|, |2j . Studies of the decays of 
quarkonia and of their decay products can also reveal im- 
portant information on QCD processes Q. The hadronic 
decays of the narrow quarkonia, states which are below 
the threshold for open flavor production, are dominated 
by couplings to gluons and the fragmentation process into 
light hadrons. The decay properties of charmonia, which 
have a relatively low multiplicity particle content, have 
been extensively studied [4]. However, little is known 
about the final state contents of bottomonia. In particu- 
lar, scarcely any experimental information exists on the 
decays of bottomonium to open charm. The CLEO Col- 
laboration has observed [j| charm production in the de- 
cays of the Xb states with branching fractions of the order 
of 10%. The ARGUS Collaboration searched @ for the 
decay T(1S) — > D* ± X and set a limit on its branching 
fraction of B < 1.9% at 90% confidence level. 



In this article, we report a study of the inclusive pro- 
cess T(1S) — > D*^X, yielding the decay branching frac- 
tion and the D*^ momentum spectrum in the T(1S) rest 
frame, using data recorded by the BABAR Collaboration 
at the Y(2S) resonance. The decay Y(1S) -> D*±X can 
proceed through the QED virtual photon annihilation 
process, T(1S) — » 7* — > cc, followed by the hadronization 
of the cc system. The expected decay rate and the D*^ 
momentum spectrum from this process can be accurately 
estimated from the measured properties of the T(1S) de- 
cays and the charm fragmentation function measured at 
the center-of-mass energy yfs ~ 10 GcV. Other QCD 
processes such as the splitting of a virtual gluon Q Hal o r 
the annihilation of the bb system in an octet state [lOj . 
have also been suggested as major contributors to this 
decay channel. Measurements of the D*^ yield and of 
its momentum spectrum can help test the predictions of 
the proposed QCD mechanisms, and possibly reveal the 
presence of new physics processes with exotic couplings 
to heavy quarks [111, EH • 
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II. THE BABAR DETECTOR 

The results presented in this work are based on data 
collected at center-of-mass energy corresponding to the 
mass of T(2S) resonance with the BABAR detector at the 
PEP-II asymmetric energy e + e _ storage ring operating 
at the SLAC National Accelerator Laboratory. The data 
consist of 14.4 fb _1 of integrated luminosity, correspond- 
ing to 98.6 ± 0.9 million T(2S) mesons produced. The 
study of T(1S) decays is performed by reconstructing the 
decay chain T(2S) — > tt + tt~Y(1S), which yields approx- 
imately 17.8 million T(IS) decays. An additional off- 
resonance data sample corresponding to 44.5 fb _ col- 
lected at y^s about 40MeV below the T(AS) resonance 
is used to study the background. A GEANT4-based 
simulation of the detector is used to determine the prop- 
erties of the signal and to study the background sources. 

A detailed description of the BABAR detector can be 
found elsewhere [14j . The tracking system is composed 
of a 5 layer silicon vertex tracker (SVT) and a 40 layer 
drift chamber (DCH) in a 1.5 T magnetic field. The 
SVT provides a precise determination of the track im- 
pact parameters and angles near the interaction point 
(IP) with 15 /im spatial resolution at normal incidence at 
a radius of 3.2 cm, and is capable of stand-alone track- 
ing for low momentum particles down to 50 MeV/c of 
transverse momentum p t . The DCH, together with the 
SVT, provides a precise measurement of the momenta 
and azimuthal angles of charged particles with a reso- 
lution (Jpjpt — (0.13 p t © 0.45)%, where p t is in units 
of GeV/c. Charged hadron identification is achieved 
through measurements of the specific ionization energy 
loss in the SVT and DCH, and of the Cherenkov angle 
from a detector of internally reflected Cherenkov light 
(DIRC). A CsI(Tl) electromagnetic calorimeter (EMC) 
provides photon detection, electron identification, and 
7T°, 77 and Kl reconstruction. Finally, the instrumented 
flux return (IFR) of the magnet allows discrimination of 
muons from pions and detection of neutral kaons. 



III. CANDIDATE RECONSTRUCTION AND 
SELECTION 

The T(2S) —> tt + it~Y(IS) candidates are identified by 
forming pairs of oppositely charged tracks whose recoil 
mass is consistent with the mass of the T(1S) resonance, 
when the tracks are interpreted as pions. The recoil mass 
-^recoil ls computed using 

Mrccoil = V(Pe+e- ~ P^f (1) 

where P e + e - is the known 4-momcntum of the e + e~ sys- 
tem and P,^ is the reconstructed 4-momentum of the 
7r + 7r~ pair. The pion tracks are required to have en- 
ergy losses and Cherenkov angles consistent with the pion 
hypothesis. The track pair is fitted to a common ver- 
tex and the probability of the vertex fit is required to 
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be greater than 1%. The measured di-pion mass dis- 
tribution peaks near 0.52GeV/c 2 for signal events [l5j . 
whereas background events are approximately uniformly 
distributed in the kinematically allowed mass interval 
[0.28,0.56] GeV/c 2 . Requiring the mass of the pion pair 
to be greater than 0.4GeV/c 2 retains 96% of the signal 
candidates while rejecting approximately 1 /3 of the back- 
ground events. Figure [T] shows the recoil mass distri- 
bution for the event sample passing the above selection 
criteria; a signal region consisting of two standard devia- 
tions around the T(1S) mass is highlighted (cross hatch- 
ing) , as well as two sideband regions used for background 
studies, the lower ([9432.1, 9444.3] MeV/c 2 ) and upper 
([9477.7, 9490.0] MeV/c 2 ) sidebands (diagonal shading). 
These events form the full event-set used in the mea- 
surement of the D*^ yield. 

We reconstruct D candidates using the decay chain 
D*+ -> D°n + , D° -> K-JT+. A kaon candidate and 
an oppositely charged pion candidate are combined to 
form the D° candidate. The identification efficiency for 
kaons (pions) is about 98% (93%); the misidentification 
rate of kaons (pions) as pions (kaons) is about 5% (15%). 
The identification performance is obtained from a con- 
trol sample of inclusive D* + -t D a -K + ,D° -> K~tt + . 
The kaon and pion tracks are geometrically constrained 
to originate from a common vertex and the probabil- 
ity of the vertex fit is required to be greater than 1%. 
The mass of the D° candidate is required to be within 
75 MeV/c 2 of the nominal D° mass, which corresponds 
to about 18 times the experimental resolution on the 
D° candidate mass. This large mass interval is neces- 
sary for the subtraction of the combinatorial background. 
The D° candidate is finally combined with a soft pion 
with its charge opposite to that of the kaon candidate 
to form a D* + candidate. The mass difference between 
the D* + and the D° (Am) is required to be in the in- 
terval [143.20, 147.64] MeV/c 2 , which corresponds to ap- 
proximately six times the experimental resolution. The 
soft pion and the D° candidates are fitted to a common 
vertex constrained to originate from the interaction re- 
gion. The probability of the D* ± vertex fit is required 
to be greater than 1%. For events with multiple candi- 
dates, the candidate with the best combined vertex fit 
X 2 , defined as the sum of the x 2 values from the vertex 
fits described above, is kept. The multiplicity of the re- 
constructed candidates in simulated signal MC events is 
1.2, after the final selection. 74% of these candidates are 
correctly matched to a signal candidate. The best can- 
didate algorithm retains 90% of the correctly matched 
candidates and 68% of the ones not correctly matched. 

IV. SIGNAL EXTRACTION 

The sample of D* ± candidates is studied in intervals 
of the scaled momentum x p , defined as: 



where pz>*± is the D momentum in the rest frame of 
the T(1S), J3 max = jE^ ax -m 2 Dt , + , E max = m r(ls) /2 

and riiD*+ is the world average of the D* + mass (l6| . 

The sample is divided into x p intervals of 0.05 width 
in the range [0.1, 1.0]; the region x p < 0.1, which is dom- 
inated by combinatorial background, is excluded. 

The invariant mass distribution of the D° candidates 
in each x p interval is used to determine the D*^ yield 
from T{IS) -> D^X. The D° mass distribution is 
obtained from the K~tt + candidates mass distribution 
by two background subtractions. Combinatorial back- 
grounds, events that are not T(2S) — > it + tt~Y(1S) de- 
cays, are removed by subtracting the lower and upper 
sidebands of the tt + tt~ recoil mass. The K~ir + invari- 
ant mass distribution from the sidebands is rescaled to 
the expected number of background events in the sig- 
nal region to determine the K~tt + mass distribution 
from the combinatorial background component under the 
T(1S) peak. In addition, the K~tt + mass distribution for 
"wrong-sign" D°(-^ K~ir + )n~ combinations (where the 
soft pion has the same charge as that of the kaon can- 
didate) is used to subtract the D* combinatoric back- 
ground including a possible peaking backgrounds from 
D°(— > K ~ir + )ir + combinations, involving a true D° de- 
cay and a random soft pion. This method leads to a small 
over-subtraction of signal events due to doubly Cabibbo 
suppressed (DCS) D° decays reconstructed as wrong-sign 
combinations. This is accounted for in the final esti- 
mation of the branching fraction. The background sub- 
tracted invariant mass distribution of D° candidates in 
the full x p range is shown in Figure [2l 

Finally the invariant mass distribution of the D° candi- 
dates in each x p interval is fitted to a probability density 
function (p.d.f.) using a minimum \ 2 estimator. The fit- 
ted p.d.f., P(m), is the sum of a signal p.d.f., P s i g (m), and 
a p.d.f. which accounts for unsubtracted backgrounds, 

-fbkgM, 

P(m) = n sig x P sig (m) + n hkg x P hkg (m) (3) 

where n s i g and ribkg are the number of signal and back- 
ground events in the fitted region. The fit region corre- 
sponds to the D° mass range in which we accept signal 
candidates [m D o ~ 75 MeV/c 2 , m D o + 75 MeV/c 2 ] . The sig- 
nal p.d.f. is the sum of two Gaussian functions with the 
same mean: 

Psi s (m;f,n, cri,<r2 ) = /G(m;/z, (Ti) + (1 - /)G(m;/i,0 2 ) 

The background p.d.f. is a linear function: 

Pbk g (m;fi,pi) = 1/w +pi(m - fx) (5) 

where w is the fit range. The parameters of the signal 
p.d.f., oi, 02 and / are determined from a fit to the cor- 
responding distribution from Monte Carlo (MC) simula- 
tion. However, the mean of the D° mass, //, is fixed to the 
value determined from a fit to the D° mass distribution 
in the full x p interval. 
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FIG. 1: Distribution of the recoil mass, M leco u, for the se- 
lected T(2S) — > ir + ir~T(lS) events. The cross hatching 
shows the signal region, and the lower and upper sideband 
regions are indicated by the diagonal shading. 



FIG. 3: Reconstruction efficiency for the decay chain 
T(2S) ->■ 7r + 7r-T(lS) ) Y(1S) -> D^X as a function of the 
scaled D* ± momentum x D . 
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FIG. 2: Distribution of the D° invariant mass for the com- 
plete [0.1;1.0] x p range after subtraction of the combinatoric 
background and wrong-sign combinations. The solid line rep- 
resents the fit to the data of the p.d.f. described in the text. 



The event selection efficiency is determined using a 
simulation study of the signal and background processes. 
Signal events are obtained by generating the transition 
T(2S) —> n + n~Y(lS) according to the decay model de- 
termined by CLEO [lH, followed by the decay r(lS)_^ 
cc and the hadronization of the cc pair via JETSET [l7| ■ 
Signal events are required to contain at least one I?* ± me- 
son after the hadronization process. The small fraction 
of events (0.4%) containing both D* + and D*~ decays is 
accounted for by normalizing the efficiency to the num- 
ber of signal decays generated. The selection efficiency as 
a function of x p , e(x p ), is shown in Figured] The depen- 
dence on x p is mainly due to the reconstruction efficiency 
of the slow pion from the D*^ decay. The average re- 



construction efficiency in data depends on the measured 
x p distribution and can be estimated from the relation 

< ^ta >= s J:;^i p) = (17.7 ±0.3)% where the 

error is statistical only. The ratio of the % 2 to the number 
of degrees of freedom for the individual fits ranges from 
0.5 to 2.5, with 16 degrees of freedom. 



V. RESULTS 

Figure |4] shows the efficiency-corrected distribution of 
the D* ± yield as a function of x p . The branching fraction 
for the inclusive decay T(1S) —> D*^X in the x p range 
[0.1, 1.0] is computed from: 



B[T{1S) -> D^X] 



fcDCS x B, 



decay 



x N- 



T(1S) 



(G) 



(2.52 ± 0.13(stat) ± 0.15(syst))% 



where n s i g — E Xp n s i g (x p )/e(x p ) = 11845 ± 596 is the 
efficiency-corrected signal yield in the x p range [0.1, 1.0], 
^dcs = (99.62 ± 0.02)% is a correction factor to account 
for the subtraction of DCS D° decays, 2?d e cay is t ne prod- 
uct of the branching fractions [ll| in the D decay chain 
B[D*+ -> D°tt + ] = (67.7±0.5)% and B[D° -> K L n+] = 
(3.91±0.05)%, N r(2S) = (98.6±0.9) x 10 6 , and N r(is) = 
N T (2S) x B[Y(2S) -> 7r+7r-T(15)] = (17.8 ± 0.4) x 10 6 
is the number of T(1S) mesons produced in this decay 
chain. 

We verify that our analysis procedure is unbiased by 
fitting off-resonance data and a Monte Carlo simulation 
of the background; we find no significant signal. We also 
compare the lower and upper sidebands and use the D° 
mass sidebands instead of the recoil mass to subtract the 
background, and we find no significant shift in the sig- 
nal. The sources of systematic uncertainty are listed in 
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FIG. 4: Signal yield as a function of x p . The solid line repre- 
sents the expected contribution from the virtual photon pro- 
cess pH ]. 



TABLE I: Summary of the systematic uncertainties on 
B[Y(1S) ->■ D^X]. 



Sources of systematic uncertainty Relative uncertainty 



Slow ir ± reconstruction efficiency 


3.0% 


Afrccoii selection 


2.8% 


Sdccay X B[T(2S) -> 7Y + 7Y-T(1S)} 


2.3% 


Generated x p distribution 


2.2% 


PID 


1.6% 


Tracking efficiency (excl. slow pion) 


1.6% 


T(2S) decay model 


1.2% 


T counting 


0.9% 


Background curvature 


0.4% 


MC efficiency 


0.4% 


Signal shape 


0.3% 


fcDCS 


0.02% 


Total systematic uncertainty 


5.9% 



Table HI The main contributions come from the uncer- 
tainties in the knowledge of the slow pion reconstruction 
efficiency and the selection efficiency of T(1S) decays in 
the recoil mass signal region. The former is determined 
from a control sample of D* + — > D°ir + decays by com- 
paring the efficiency in data with that in MC events, 
for the soft pion momentum range [50,400] MeV. The 
efficiency is extracted from a study of the angular dis- 
tribution of the soft pion in the rest frame of the D* 
meson. The M reco ii selection systematic uncertainty is 
obtained by comparing the recoil mass distribution for 
signal events in the full x p range [0.1,1.0] in data, with 
the distribution in Monte Carlo simulated events. The 
fit to data with the sum of two Gauss functions gives 
an r.m.s. of 2.9 MeV while the fit to MC events gives 
3.3 MeV. The efficiency is estimated from the integral of 
the fitted function in a window around the T(lS) mass 
of ±2x the r.m.s. on MC (the recoil mass signal region). 
The efficiency in data is 96.3% while in MC events is 



93.6%, which corresponds to a relative systematic error 
on the result of 2.8%. The uncertainty associated with 
the generated x p distribution is determined by reweight- 
ing simulated signal MC events according to the x p distri- 
bution measured using data. In addition the parameters 
of the T(2S) decay model have been varied within their 
uncertainty and the resulting relative efficiency variation 
has been taken as the systematic uncertainty. The un- 
certainty in the particle identification efficiency (PID) is 
derived from a study of a <fr — > K + K~ control sample 
and by removing the PID requirement from the selection. 
The dominant systematic uncertainties in the T count- 
ing come from the modeling of the track reconstruction 
efficiency and of the total energy of the events. The sig- 
nal shape uncertainty is due to data-MC differences in 
the D° mass signal distribution. A possible curvature of 
the background is extracted from off-resonance data, and 
the systematic uncertainty is obtained by adding the cor- 
responding second order polynomial to the background 
p.d.f.. The uncertainties due to MC efficiency, fees an d 
-Bdecay x B[Y(2S) — > 7r + 7r~ Y(1S)} arise from imperfect 
knowledge of these parameters. 

VI. DISCUSSION AND CONCLUSION 

Figure @] shows the expected x p distribution for Z3* ± 
production from the QED virtual photon annihilation 
process, T(1S) — > 7* — > cc. The shape is obtained 
from the measured D*^ fragmentation function at y/s — 
10.5 GeV [18| and the normalization is computed from: 

B[T(1S) -> 7* -J- D*±X] x R had (7) 

(Jqq 

x B[Y{IS) ->■ 

where, i?had = cr(e + e~ — > hadrons)/er(e + e~ — > = 
3.46 ± 0.130, B[T(1S) = (2.48 ± 0.05)%[l6j], 

and = (17.7 ± 2.2)%[l6j is the measured D* ± yield 

from e+V ->• qq at ^/s = 10.5 GeV. We find B[T{1S) -> 
7* -> D*±X] = (1.52 ± 0.20)% . 

Our measured branching fraction exceeds the expected 
rate from the QED virtual photon process from Eq. [7] 
by (1.00 ±0.28)% (including the systematic uncertainty) 
which corresponds to 3.6 standard deviations. While 
the measured x p spectrum agrees in shape with that 
of the virtual photon process for x p > 0.75, there is a 
significant excess for x p < 0.75. The probability that 
the measured spectrum is consistent with the expected 
distribution from the virtual photon, normalized using 
Eq. ([7]), is 1.2 x 10~ 5 confidence estimated from a binned 
X 2 test. The excess is compatible with the contribu- 
tion expected Q from the splitting of a virtual gluon, 
(1.20 ± 0.29)%. This does not leave much room for the 
octet contribution [l(| , which is also disfavored from the 
shape of the excess as a function of x p . 

In summary, using the data collected with the BABAR 
detector at the T(2S) resonance, we have observed for 
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the first time the decay of T(1S) mesons to open charm. 
We have measured the branching fraction B[T(1S) — > 
D*±X] = (2.52 ± 0.13(stat) ± 0.15(syst))% and the D*± 
momentum distribution in the rest frame of the Y(IS). 
We find evidence for a significant excess of D*^ produc- 
tion with respect to the expectation from the virtual pho- 
ton annihilation process. 
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